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> : ABSTRACT 
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^: 

^ I We have simulated cosmic microwave background (CMB) anisotropy maps 

I I for several COi? -E-DMR-normalized cold dark matter (CDM) cosmogonies, 

J> ' to make predictions for the MAP experiment, an upcoming whole-sky CMB 

i/^ [ anisotropy space mission. We have studied the sensitivity of the simulated MAP 

O I data to cosmology, sky coverage, and instrumental noise. Given an accurate 

knowledge of instrumental noise, MAP data will ably discriminate amongst 
the cosmogonies considered, and superbly determine the topology of the initial 
fluctuations. 

Ph| a correlation function analysis of the simulated MAP data results in a 

O ' very accurate measurement of the acoustic Hubble radius at decoupling. A 

I low-density open CDM model with fio = 0.4 can be distinguished from the 

' Qq = 1 fiducial CDM model or a spatially-fiat CDM model with a cosmological 

' constant and = 0.4 with more than 99% confidence from the location of the 

^ ' acoustic "valley" in the correlation function. 

c3 I A genus analysis of the simulated MAP data indicates that in cosmogonies 

with Gaussian random-phase initial conditions, the horizontal shift of the 
zero-crossing point of the genus curve near the mean temperature threshold 
level z/ = should not exceed |Az/| = 0.01 (0.04) when the total effective FWHM 
smoothing is 0.3° (1.0°). The asymmetry of the genus curve at the positive and 
negative threshold levels should not exceed lAg/g^u = ±1)| = 0.8% (4%) at 
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0.3° (1°) FWHM smoothing. Deviations of the observed MAP data in excess of 
these small values will be evidence for non-Gaussian behavior. 

The amplitude of the genus curve is a measure of the shape of the 
power spectrum at the smoothing scale. Even with the expected amount of 
instrumental noise and partial sky coverage (due to the Galaxy), the MAP data 
should allow discrimination amongst the cosmogonies considered at more than 
99% confidence solely from a genus amplitude analysis. 

Subject headings: cosmology — cosmic microwave background: anisotropy — 
large-scale structure of the universe 

1. Introduction 

Following the COBE-DMK detection of anisotropy in the cosmic microwave 
background on large angular scales (Smoot et al. 1992; Wright et al. 1992; Bennett et 
al. 1996; Gorski et al. 1996), there have been many measurements of CMB anisotropy 
on angular scales down to ~ 10' (Ganga et al. 1994; Gutierrez et al. 1997; Piccirillo 
et al. 1997; Netterfield et al. 1997; Gundersen et al. 1995; Tucker et al. 1997; Piatt 
et al. 1997; Masi et al. 1996; Lim et al. 1996; Cheng et al. 1997; Griffin et al. 1997; 
Scott et al. 1996; Leitch et al. 1997; Church et al. 1997, see Page 1997 for a review). In 
addition to the DMR experiment, CMB anisotropy maps have been produced from four 
other experiments (Ganga et al. 1994; White & Bunn 1995; Scott et al. 1996; Tegmark 
et al. 1997). CMB anisotropy observations are beginning to test cosmogonical models and 
provide interesting constraints on cosmological parameters in these models (e.g., Bunn 
& Sugiyama 1995; Gorski et al. 1995, 1998; Stompor, Gorski, & Banday 1995; Ganga, 
Ratra, & Sugiyama 1996; Ganga et al. 1997a,b, 1998; Hancock et al. 1997; Bond & Jaffe 
1997; Stompor 1997; Lineweaver & Barbosa 1997; Ratra et al. 1997b). However, definitive 
results, including significant constraints on cosmological parameters, must await new CMB 
anisotropy data, covering a large fraction of the sky (to minimize sample variance) and 
acquired at a variety of frequencies (to allow for control over possible non-CMB anisotropy 
foreground contamination). Of the currently available data sets, the DMR maps best meet 
these criteria. 

Plans are afoot to launch two second-generation CMB anisotropy satellite missions — 
the Microwave Anisotropy Probe (MAP)| is expected to have data in less than four years. 
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and the Planck Surveyor[] is expected to have data in less than a decade. These missions 
are expected to measure the CMB anisotropy on angular scales larger than a fraction of 
a degree, at a variety of frequencies. MAP will probe the anisotropy at frequencies of 22, 
30, 40, 60, and 90 GHz at expected angular resolutions of 0.93, 0.68, 0.47, 0.35, and 0.21 
degrees (FWHM beamwidths), respectively. At 90 GHz the MAP sensitivity is expected 
to be ~ 35 /xK per 0.3° x 0.3° pixel. Planck Surveyor is expected to measure the CMB 
anisotropy at frequencies ranging from 30 to 900 GHz with angular resolution ranging from 
0.5 to 0.075 degrees. 

It is now possible to make accurate predictions for the CMB anisotropy, as a function of 
cosmological parameters, in many different cosmological models. In conjunction with such 
theoretically computed model CMB anisotropy spectra, the data from these experiments 
should lead to an accurate determination of cosmological parameters, such as the mass 
density parameter Qq, the baryonic-mass density parameter Qb, the cosmological constant 
A, the Hubble parameter Ho, the matter power spectrum, and the optical depth to the 
last scattering surface (e.g., Spergel 1994; Jungman et al. 1996; Bond, Efstathiou, & 
Tegmark 1997; Zaldarriaga, Spergel, & Seljak 1997; Dodelson, Kinney, & Kolb 1997). At 
the least, the data from these space missions should in principle be able to rule out most 
cosmological models currently under consideration. The present situation in cosmology 
— where one is able to make definite theoretical predictions which can be judged by 
near-future observations — is quite unusual. 

Prior to the COBE-DMR measurement of the large-scale CMB anisotropy, Gott et al. 
(1990) had generated simulated all-sky CMB anisotropy maps predicted in the spatially-flat 
r^o = 1 fiducial CDM model with Gaussian, adiabatic, scale-invariant matter fluctuations. 
The DMR data is indeed consistent with such a CMB anisotropy spectrum (Gorski et al. 
1996, 1998). 

In this paper we generate simulated all-sky CMB anisotropy maps appropriate for MAP 
in several representative cosmological models normalized to the DMR data. Our predictions 
can be directly compared to the MAP data when it is available. Hinshaw, Bennett, & Kogut 
(1995) have generated all-sky CMB anisotropy maps expected in an Qq = 1 CDM model, at 
angular resolutions of 0.5° and 1°, to study the effects of sky coverage, angular resolution, 
instrumental noise, and unresolved features in the data. Novikov & J0rgensen (1996) have 
simulated 10° x 10° sky anisotropy maps of = 1 CDM models with Qb = 0.1 and 0.03 at 
0.5°-l° resolutions. They have examined various statistical properties of the temperature 
fluctuation field, such as clustering of peaks and percolation. We focus here on generating 
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mock CMB anisotropy maps appropriate for the specific experimental configuration of 
MAP. We also use these CMB anisotropy maps to examine the ability of the expected MAP 
data to discriminate between different cosmological models. In what follows we assume that 
the MAP data will be limited by statistical noise rather than foregrounds or systematic 
errors. 

In Section 2 we present the cosmological models we study, and explain how the 
simulated CMB anisotropy maps are generated. In Section 3 we compute the CMB 
anisotropy correlation function for each of the simulated maps, and show that it is sensitive 
to the acoustic peaks in the CMB anisotropy angular power spectrum. In Section 4 we study 
the genus statistic of the generated maps, and explore its ability to discriminate between 
cosmological models and to detect non-Gaussian fluctuations in the CMB anisotropy field. 
The effect of instrumental noise is considered. Conclusions follow in Section 5. 



2. Simulation of CMB Anisotropy Maps 

To generate simulated CMB anisotropy maps for MAP we have chosen four 
representative Gaussian, adiabatic cosmogonical models, with distinct CMB anisotropy 
angular power spectra: (1) The fiducial CDM model (FCDM) with D,q — 1, the Hubble 
parameter Hq divided by 100 km/sec/Mpc h — 0.5, and with no reionization. (2) A 
spatially-fiat, low-density CDM model (ACDM) with a cosmological constant A, = 0.4, 
and h = 0.6. (3) A spatially-open low-density CDM model (OCDM) with no A, = 0.4, 
and h = 0.65. (4) A reionized fiducial CDM model (RCDM) with = 1, h = 0.5, and 
fully reionized at redshift z = 100 (i.e., with electron ionization fraction Xe = 1 aX z < 100 
and so optical depth r = 1). For all models, the baryonic-mass density parameter is 
Qb = 0.0125/i~^. Further details are given in Ratra et al. (1997a) and Sugiyama (1995). 

The FCDM and ACDM models considered here assume a scale-invariant primordial 
energy-density power spectrum (Harrison 1970; Peebles & Yu 1970; Zel'dovich 1972), as is 
predicted in the simplest spatially-flat inflation models (Guth 1981, also see Kazanas 1980; 
Sato 1981). The FCDM model is inconsistent with galaxy clustering data (e.g., Maddox et 
al. 1990; Saunders et al. 1991; Vogeley et al. 1992; Park et al. 1994), but stiU serves as a 
useful fiducial for cosmological model comparison. ACDM models (Peebles 1984) became 
popular when it was discovered that the large-scale galaxy distribution was of the form 
expected in CDM models with ^oh ~ 0.2 - 0.3 (e.g.. Park 1990, 1991; Vogeley et al. 1992; 
da Costa et al. 1994; Peacock & Dodds 1994), since these models could be accommodated 
in the then-dominant spatially-flat inflation scenario. However, an open inflating universe 
can also be created as a single bubble during the decay of an initial metastable inflation 
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state (Gott 1982, 1986). The OCDM model considered here assumes the primordial 

energy-density power spectrum of the simplest open-bubble inflation models (Ratra & 
Peebles 1994, 1995; Bucher, Goldhaber, & Turok 1995; Yamamoto, Sasaki, & Tanaka 1995). 
The low-density models, OCDM and ACDM, are consistent with most current large-scale 
structure and CMB anisotropy observations (e.g., Gott 1997; Turner 1997; Cole et al. 1997; 
Gardner et al. 1997; Croft et al. 1997; Coles et al. 1997; Jenkins et al. 1997). 

A CMB temperature fluctuation fleld on the surface of a sphere (i.e., the sky) can be 
decomposed into a sum of spherical harmonics, 6T/T{6,(f)) = a™^'"(^) 0)) where {9,(f)) 
are the angular coordinates on the sky. If the initial matter fluctuation is a Gaussian 
random field, as predicted in simple infiation models, the expansion coefficients a™ are 
independent Gaussian random variables with zero means and variances = Ce, where 

Ce is the CMB anisotropy angular power spectrum (see below for more details). 

The CMB anisotropy angular power spectra of the four models considered here are 
shown in Figure 1 (Ratra et al. 1997a; Sugiyama 1995). Following Gorski et al. (1998) and 
Stompor (1997), the angular power spectra are normalized to the arithmetic mean of the 
±2-(7 COBE-DMR 53 and 90 GHz data from the two extreme data sets: (1) galactic-frame 
maps accounting for the high-latitude Galactic emission correction and including the i — 2 
moment in the analysis; and (2) ecliptic-frame maps ignoring the high-latitude Galactic 
emission correction and excluding the £ — 2 moment from the analysis. 

At large angular scales the CMB anisotropy angular power spectra in the four models 
are of approximately the same amplitude. At smaller angular scales acoustic oscillations 
due to the pressure of the photon-baryon fluid and matter velocity perturbations at photon 
decoupling cause the angular power spectra to rise. At still smaller angular scales photon 
diffusion and the flnite thickness of the last-scattering surface smears out the fluctuations. 
These effects give rise to the characteristic acoustic peaks in the angular power spectra on 
intermediate angular scales. In an open model the geometrical effect that the circumference 
of a circle is larger than 27r times its radius moves these peaks to smaller angular scales 
relative to the spatially-flat case (Kamionkowski, Spergel, & Sugiyama 1994). In the 
spatially-flat — 1 FCDM and low-density ACDM models, the first acoustic peak is 
at £ = 217 and £ — 222, respectively. The amplitudes of the peaks in these models are 
different. In the open case the first peak is at £ = 349, while in the reionized RCDM model 
the peaks are smeared out. Consequently, the presence, location, and amplitude of the 
acoustic peaks are a valuable probe of the geometry and content of the universe. 

In real observations and analyses, the temperature fluctuations on the sky are convolved 
with the beam pattern of the experiment, and are often further smoothed to increase the 
signal-to-noise ratio. If the beam pattern is B{9) and the smoothing filter is F{9), the 
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observed CMB anisotropy temperature field is 

oo t 

= E E «r>^r(^, miFt. (i) 

where and are the coefficients of the beam pattern and the smoothing filter expanded 
in Legendre polynomials, and wc consider only harmonics higher than the dipole. Assuming 
the idealized case of an all-sky data map, the inverse relation for the harmonic amplitude is 

af = J dn^-^{e, <i^r/BeFi. (2) 

For a Gaussian beam (or a smoothing filter) with a width A <^ 1 to a very good 
approximation, 

S^ = exp(-(£+l/2)2AV2). (3) 

We use equation (1) to simulate the CMB anisotropy maps for MAP. We limit the 
summation to £ < 1500. For each cosmogony the real and imaginary components of a™ are 
randomly sampled Gaussian variates with zero means, and variances equal to Q/2 for the 
m 7^ harmonics and equal to Q for the m — harmonics. We adopt a Gaussian beam 
with FWHM 9b = 0.21° (or width A = 9b /^V^ In 2 in equation [3]). The temperature 
fiuctuations are then computed at 1801 x 3600 points in {9, 0) space on a grid with 0.1° 
spacing. 

The data from MAP is expected to include instrumental noise of ~ 35 /jK per 
0.3° X 0.3° pixel. If the noise is white, its power spectrum is constant Ce = C^. Consider a 
pure noise map convolved with a top hat filter F with area of 0.09 square degrees. Then the 
variance of the noise temperature fluctuation, which is set to (35 fiK/To)'^ per 0.09 square 
degrees, is given by equation (5) below (at = with replaced by Cn)- The top hat 
filter is F{6) = l/TiO'jj for 6 < 6h, and otherwise. The coefficients of its expansion into 
Legendre polynomials are, for 6'^^ -C 1, 

oo / \n 

Fi = T. ns^rrvdi^ + imonm'-. (4) 

^0 {niy{n+l) 

We obtain from equation (5) Cn — 5 x 10"^^ for MAP. One can study the effects of noise 
on statistics using this noise power spectrum. In practice, of course, the MAP results will 
include an estimate of the noise at each pixel derived from the many observations taken 

as the satellite's horns sweep out complicated paths on the sky. The satellite horns are 
paired, so subtracting their signals gives the noise map; this allows an extremely accurate 
characterization of the noise as a function of position on the sky. Such a detailed noise map 



- 7- 



will allow a much more careful accounting of instrumental noise. (For example, the noise at 
the ecliptic poles is expected to be about half of that at the echptic plane) . 

The maps resulting from the ACDM and OCDM models with simulated (constant 
over the sky) instrumental noise, projected according to the MoUweide projection method 
(Snyder 1993), are shown in Figure 2 (Plate 1). These maps are smoothed with a 0.214° 
FWHM Gaussian, which makes the total smoothing scale of the signal 0.3°. The Mollwcide 
equal area projection shows the whole sky in a single map, but structures are distorted. 
In the topology study we use conformal stereographic projection maps where shapes of 
structures are preserved locally. Stereographic maps of the CMB fluctuations simulated 
for the OCDM and ACDM models and convolved with a 0.978° FWHM Gaussian (which 
makes the total smoothing scale of the signal 1.0°) are shown in Figure 3 (Plate 2). At this 
smoothing length most structures in the map generated by MAP will correspond to actual 
CMB anisotropy fluctuations (see Figure 8). It is to be noted that the OCDM and ACDM 
maps are visually very similar — this is a consequence of the expected instrumental noise 
and smoothing (see Figures lb and Ic). We emphasize, however, that this does not mean 
that MAP can not statistically distinguish between these models. 

3. Correlation Function 

Differences between cosmological models are often displayed in terms of the angular 
power spectrum, as in Figure 1. Most models have relatively small differences at large scales 
{£ < 30), but they often have relatively large differences at small scales {£ > 100). The 
structure of the acoustic peaks at smaller angular scales depends on the cosmogony and the 
cosmological parameters. Hence, to test cosmogonies and measure these parameters, it is 
essential to have an accurate measurement of the location and amplitude of the acoustic 
peaks. Note that the acoustic peaks cause an oscillatory wave in the power spectrum (see 
Figure 1). Furthermore, the amplitudes of individual harmonics in the power spectrum 
are destined to fluctuate with respect to the mean. To measure precisely the amplitude 
and location of the acoustic peaks in the power spectrum one must make a fit over many 
harmonics. 

An important integral of the power spectrum is the two-point auto-correlation function, 
the real-space Fourier counterpart of the power spectrum. Since the hot and cold spots due 
to the acoustic peaks in the power spectrum are compact in real space (their size is fixed 
by the acoustic Hubble radius at decoupling — see Figures 2 and 3), it is useful to explore 
them in real space. (Note that in Fourier space the peaks are well- approximated by a wavy 
sine- function.) While the correlation function and the power spectrum are simply related. 
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it is useful to analyze data both in real space and in Fourier space (e.g., Hinshaw et al. 
1996; Gorski et al. 1996; Wright et al. 1996). 

The correlation function of CMB temperature fluctuations, corresponding to equation 
(1), is 

C(^) ^ ^Y^Cei2e+l)Piicose)BjF^. (5) 

Here are Legendre polynomials and we have ignored the monopole and dipole harmonics. 
The power spectra shown in Figure la can be approximated by an exponential minus 
a sinc-function. For example, the acoustic peaks of the FCDM power spectrum can be 
approximately described by the function —smc{4.5£ /imax)- Consequently, these acoustic 
peaks should appear as a negative top hat in the correlation function. 

Figure 4 shows the correlation functions of the four models considered. Thick curves 
are the correlation functions transformed from the power spectrum curves in Figure 1. Thin 
solid curves, tracing the thick ones at large angular separations, account for the convolution 
with a 0.21° Gaussian beam in Figure 4a and with additional 0.214° Gaussian smoothing 
in Figure 46. The acoustic peaks in the power spectrum indeed appear as a valley in this 
plot. The position, depth, and width of the acoustic valley depend on the model and the 
cosmological parameters. We have generated seven sets of a^^'s for the FCDM model, and 
have computed their correlation functions. These are also plotted, using dots, in Figure 4. 
At separations > 1° cosmic variance causes a large scatter amongst the seven correlation 
functions. Hence it is very difficult to discriminate between models using large angular 
scale data alone. However, direct measurement from these realizations shows that the 
FCDM acoustic valley in the function e^/'^C{e) is consistently located at ^ = 1.13° ± 0.02°, 
comparing excellently with the theoretically predicted FCDM acoustic valley location at 
e = 1.13° for a 0.21° beam (and at 1.17° in the unsmoothed map). The OCDM model 
acoustic valley lies at a smaller separation of 0.65° ± 0.03° (0.69° in the unsmoothed map) 
— clearly distinguishable from the acoustic valleys of the FCDM and ACDM models (some 
13-(7 difference for the fixed cosmological parameter values used here) . 

On the other hand, the spatially-fiat FCDM and ACDM models have acoustic peaks at 
nearly the same angular scales. Therefore, they have to be differentiated by the amplitude 
of the power spectrum or the correlation function. It can be noted from Figure 4 that the 
correlation functions of the FCDM and ACDM models are statistically significantly different 
at separations fa 0.5°. In the case of an ideal sample with whole-sky coverage and with no 
noise, the MAP data would distinguish the correlation function of the FCDM universe from 
that of the ACDM case at the 2.6-cr level at angular separations near 9 = 0.5°. Note that 
this ignores the relative uncertainty between the DMR normalizations of the FCDM and 
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ACDM models. The numerical value of this relative uncertainty is not known, however, the 
absolute model-dependent uncertainty in the DMR normalization is ~ 10-12% (in Qrms-PS) 
the quadrupole amplitude of the CMB anisotropy, Gorski et al. 1998). Also note that, for 
this comparison between FCDM and ACDM, cosmological parameters such as h and Qb 
are held fixed. 

To estimate the effects of instrumental noise on the correlation function analysis of the 
MAP data wc use the noise power spectrum, = 5 x 10~^^, derived above. This noise 
power is added to the CMB anisotropy signal, and the map is then further smoothed with a 

0. 214° FWHM Gaussian. That is, we generate a set of (a^fi^ + n)Fe up to £ = 2500 where 
n is a Gaussian random variable with zero mean and (n^) = Cjv, and and represent 
the 0.21° beam and 0.214° smoothing filter, respectively. The corresponding correlation 
function is then computed. Prom seven realizations of the noise-added correlation functions 
for FCDM and OCDM the positions of the acoustic valleys are measured to be at 

9 = 1.13° ± 0.03° and 0.65° ± 0.04°, respectively. The correlation function acoustic valleys 
of the FCDM and OCDM models are now 9-0" apart. It is also noted (from Figure 4) 
that the statistical uncertainty in the amplitude of correlation function at separations 
6 > 0.5° is hardly affected by the instrumental noise expected for MAP. This means that 
cosmic variance, rather than instrumental noise, dominates the statistical uncertainty in the 
correlation function at these scales. (This is also true for the genus statistic — see Table 

1. ) Therefore, the above argument on discrimination between the spatially-fiat models on 
the basis of the amplitude of the correlation function remains essentially unchanged. 



4. Topology 

The power spectrum and correlation function measured from the observed CMB 
anisotropy maps provides important information on the nature of density fiuctuations in 
the early universe, thus constraining structure formation mechanisms. Another important 
statistical characteristic of a CMB anisotropy map is whether or not its temperature 
fiuctuation field is Gaussian random-phase. One can test for Gaussianity of the temperature 
distribution by constructing histograms of the pixel values. To test the Gaussianity and 
phase properties of the maps the genus topology test is useful. In this section we study the 
topology of the CMB anisotropy isotemperature contour surfaces. This provides information 
on the global structure of the density fluctuation fleld. 
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4.1. The Genus 

The genus has been used as a measure of topology of the matter/galaxy distribution. 
For three-dimensional topology see Gott, Melott, & Dickinson (1986), Gott, Weinberg, & 
Melott (1987), Vogeley et al. (1994), and Cole et al. (1997); for two-dimensional topology 
see Melott et al. (1989), Park et al. (1992), and CoUey (1997). The genus has also been 
used as a measure of the topology of CMB anisotropy fields (Coles 1988; Park & Gott 1988; 
Gott et al. 1990; CoUey, Gott, & Park 1996; Kogut et al. 1996; Schmalzing & Gorski 1997; 
Winitzki & Kosowsky 1997). For the two-dimensional CMB anisotropy temperature field, 
the genus is the number of hot spots minus the number of cold spots, i.e., at a temperature 
threshold level v the genus is 

where k is the signed curvature of the isotemperature contours C. The genus curve as 
a function of the temperature threshold level has a characteristic shape for a Gaussian 
random-phase field. Changing the shape of the CMB anisotropy power spectrum only 
changes the amplitude of the genus curve. Hence the genus curve provides a good measure 
of both non-random phases and the shape of the power spectrum. The above real-space 
definition of the genus allows it to be directly used for studying the statistics of the 
acoustic peak hot and cold spots. Previous measurements of the CMB anisotropy genus 
(see references above) were consistent with the Gaussian random-phase hypothesis. 

The theoretical genus per steradian of a two dimensional Gaussian field with correlation 
function 0(6) is (Gott et al. 1990) 

where (7*^"^ = {—)'^^'^ {d^C / d9"')g=o, and the threshold temperature fluctuation is 
vVC^ — va. For a CMB anisotropy temperature map convolved with beam B and 
smoothing fllter F, C*^°^ and (7*^^^ are given by 

C^"^ = ^T.i'^^ + l)CtBlF^ (8) 
C(^) = ^Y.(^{(i+im+l)C,BjFl (9) 

where we have used the relation P^'(l) = l{l -\- l)/2 (see equation [5]). The genus per 
steradian is then 

„(.) = 1 T.t(t+mt+mBin 
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We note again that the genus curve shape is fixed by the Gaussian random-phase nature 
of the field, and its amphtude depends only on the shape of the CMB anisotropy power 
spectrum and not on its amplitude. 

For scale-free energy- density perturbations with matter power spectrum 
P{k) = (l^fcP) = in a spatially-fiat universe, the CMB anisotropy angular 
power spectrum is (Peebles 1982; Bond & Efstathiou 1987) 

^ 2-' A r(3-n)r(g^) 

We can estimate the slope of the CMB anisotropy power spectrum near the smoothing scale 
from the amplitude of the genus curve by using equation (10) and this formula. 



4.2. Topology of the MAP Simulations 

When the MAP data becomes available one can measure the genus of the observed 
CMB anisotropy and so constrain cosmological models and structure formation mechanisms. 
Here we use our simulated CMB anisotropy maps to study the sensitivity of the genus 
statistic to models in the presence of instrumental noise and when the sky coverage is 
incomplete. 

We use the stereographic projection to map the CMB anisotropy sky simulations onto 
a plane. The stereographic projection is conformal and preserves local geometry. For a 
hemisphere the projection is defined by 

p. 2 tan (^^/H^); ^ = (12) 

where p is the radius and is the position angle in the projected plane, b is galactic latitude 
and i is galactic longitude (and should not be confused with the Legendre polynomial 
index). Figure 3 (Plate 2) shows such projections for the OCDM and ACDM simulations. 

When computing the genus we exclude regions with \b\ < 30° to excise the Galactic 
plane region where the Galactic signal is likely to be important.^ We use the method of 
Gott et al. (1990) to compute the topology of the proper planar projection of the simulated 
CMB anisotropy sky. In Figure 5 points with error bars show the two-dimensional genus 



^ When the MAP data becomes available it will be possible to use a less restrictive, but more complex, 
Galactic-plane cut, as was done for the DMR data (Banday et al. 1997). 
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topology, at 1.0° FWHM smoothing, derived from the simulated MAP data (accounting 
for noise). The error bars are formal errors in the mean derived from the variance among 
the genus measurements in each of the eight octants of the sphere. The curves show the 
function g = Nh'e~'^^^'^, the form expected for a random-phase Gaussian field, where N 
has been adjusted to fit the points optimally via the least-squares method. For goodness 
of fit, the points along the curves should be compared with Student's t-variates, rather 
than Gaussian variates, since the error bars have been estimated from the data themselves 
(rather than independently). Following Gott et al. (1990) and CoUey (1997), we use the 
variable 

to estimate the goodness of fit. The sum runs over the 21 u values at which the genus gi 
is computed; (ji is the mean genus among the eight octants of the sphere and Sj^ is the 
variance. We find for the OCDM case = 23.6, and for FCDM = 26.3. Since 
does not follow a distribution, due to the non-Gaussian distributions of the gi^s, we 
generate 10^ Monte Carlo runs of x^ variables with 20 degrees of freedom, and 8 samplings 
within each degree of freedom, exactly as we expect in equation (|13D (since we have fit one 
parameter, the amplitude). In each run, we generate 8 independent (// = 0, a = 1) Gaussian 
variates within 20 independent groups. From the eight samples in each group i we estimate 
the mean [cji) and the standard deviation in the mean Sij \Jn — 1; (sj/ \Jn — 1) is thus a 
(yU = 0, cr = 1) t-variate. In the 10"^ Monte Carlo simulations, we compute as in equation 
(p!3|), by summing the squares of the twenty t- variates. We find that the values of 23.6 
and 26.3 fall at the 43% and 54% levels of the cumulative distribution of the 10"^ simulated 
X^ variables, demonstrating an excellent agreement with the theoretical curve. 

In addition to noise, several non-CMB anisotropy foregrounds could contaminate the 
MAP data. Perhaps the foreground of most consequence for the genus analysis is the radio 
point-source background (mainly AGNs). To assess its significance, we have added point 
sources to the 0.3° FWHM smoothed OCDM map. These point sources (also smoothed) 
are taken to have a number distribution in flux, /, given by dN{f) = kf~'^'^df . The 
normalization k is chosen to reflect the real sky and has an average of one source at 3 /iK or 
brighter per 0.3° x 0.3° pixel (Holdaway, Owen, & Rupen 1994). Comparing to the OCDM 
map with noise alone, we find that the point sources introduce an rms scatter in g{v) of 23 
(relative to the point-source- free map), significantly smaller than the typical uncertainty in 
(7(1/), which is 42 in the map without point sources. We have also compared the genus values 
(between the maps with and without point sources) using the x^ statistic. We find x^ = 2.8 
out of an expected 21 (21 z/ values), confirming that the effect of the point sources on the 
genus is very small. Since larger smoothing lengths would diminish the signal from point 
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sources even further (as with the noise), we conclude that the effects of the point-source 
foreground on the genus are small. 

Galactic foregrounds, including warm dust and free-free emission, have been selected 
against in choosing the frequency coverage of MAP. At the optimum frequencies, these 
foregrounds are expected to be nearly an order of magnitude smaller than the CMB 
anisotropy. Furthermore, in the analysis of the DMR data (Bennett et al. 1996; Gorski 
et al. 1996, and references therein) the COBE-DIKBE 140 /xm sky map (Reach et al. 
1995) was used to correct for dust and free-free emission contamination, and the 408 MHz 
all-sky radio survey (Haslam et al. 1982) was used to correct for synchrotron emission 
contamination. A similar approach, strengthened by new all-sky survey data, should 
(hopefully) suffice for deahng with these contaminants in the MAP data. 

Another possible foreground is the Sunyaev-Zel'dovich effect in clusters of galaxies 
(Sunyaev & Zel'dovich 1972). Here hot electrons Compton scatter the CMB photons to 
uniformly higher frequency, introducing an apparent decrease in brightness temperature 
as great as ~ 1 mK on the few arc-minute scale (for very massive clusters with a velocity 
dispersion of ~ 1000 km/s, Markevitch et al. 1992). At smoothing lengths of 0.3° and 
higher, the effect is reduced to ~ 10 //K, significantly below the peaks in the CMB 
anisotropy (see Figure 2). Also, X-ray surveys of hot gas in the intracluster medium 
can be used to locate the massive clusters causing the shift, and to help estimate the 
expected magnitude of the shift (e.g., Klein et al. 1991). (In the MAP data, the large 
Sunyaev-Zel'dovich signal will come from the Coma cluster; this region can be excised from 
the analysis. This signal can also be estimated by projecting the thermal fluctuations out 
of the two highest frequency MAP channels. A combination of these cuts should be able to 
remove most of the Sunyaev-Zel'dovich contamination from the MAP data.) 

Non-Gaussian features in the CMB anisotropy will affect the genus curve in many 
different ways. Here we consider two cases that can be quantified. Non-Gaussianity can 
shift the observed genus curve to the left (towards negative thresholds) or right, near the 
mean threshold level. It can also alter the amplitudes of the genus curve at positive and 
negative levels differently, causing \g{iy — 1)| 7^ \g{^ — ~1)|- The direction and degree of 
shift and asymmetry of the genus curve depend on the number, size and height/depth of 
hot and cold spots. We have computed limits on these non-Gaussianity measures from 
our simulated maps. In the FCDM model we find that the MAP data will detect a shift 
of the genus curve near z/ = (the five u values in Figure 5 centered on u = 0) for 
|Ai/| > 0.01 (0.04) for FWHM smoothing of 0.3° (1.0°); the limits are the 2-a error bars 
derived from the variances between fits for Au from the eight octants. The asymmetry of 
the genus amplitude near i/ — ±1 (the five i/ values centered on 1 and —1) can be detected 
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for |A^/^(i/ = ±1)1 > 0.8% (4%) at 0.3° (1°) FWHM smoothing (also 2-a). The hmits for 
\Au\ and \Ag/g(u = ±1)| are very similar for OCDM and ACDM. 

4.3. Genus Amplitude as a Check of Cosmological Model 

Prom Figure 5, it is clear that the amplitude of the genus curves for FCDM and 
OCDM differ significantly. We first consider the genus amphtude for the case when 
instrumental noise is ignored. In the four cosmogonies considered, and in the scale-invariant 
power-law model, wc expect from equation (10) a theoretical genus per stcradian of 
g{iy = 1) = 101.0 ± 6.2 (FCDM), 115.4 ± 6.9 (ACDM), 74.3 ± 6.6 (OCDM), 45.2 ± 3.9 
(RCDM), and 44.4 ±3.9 {n = 1 scale-invariant power-law model), when the total smoothing 
scale (by the beam and the filter) is 1° FWHM. The uncertainty limits here account for 
cosmic variance, computed from actual realizations of a^'s. Note that these hmits are 
computed from the Galactic-plane-excluded maps. Also, since the genus amplitude does 
not depend on the amplitude of the CMB anisotropy angular power spectrum, these results 
are insensitive to the uncertainty in the DMR normalization of the models. 

We have also computed the genus directly from simulated maps that include receiver 
noise, for a single reahzation of each model. (Since all the model computations use the 
same random number seed, the cosmic variance effect is minimal.) In Figure 6 we have 
plotted these genus amplitudes at i/ = 1 for the FCDM, OCDM and ACDM models. For 
comparison, we have also plotted the genus for the same cosmogonical simulations but now 
ignoring noise. In the inset, we have zoomed in on the values at 0.3°, 0.5° and 1.0° total 
FWHM smoothing, with 1-a error bars derived from comparison of the fits to the genus 
values within each of the eight octants (these error bars agree very well with the Ax^ = 1 
level, as expected, c.f. Press et al. 1992). We have listed in Table 1 the genus values, with 
statistical uncertainty and expected cosmic variance, from Figure 6. Note that the 0.3° 
column values are subject to finite pixehzation of both the model MAP results and the 
projected maps (Melott et al. 1989). Tests have shown that amphtude ratios between 
models remain very nearly constant when moving to higher resolution in the projected 
maps. Hence, it is the finite pixehzation of our simulations which ultimately limits the 
resolution for the genus computation. 

The three cosmogonies of Figure 6 and Table 1 have significantly different genus 
amplitudes, even when the expected instrumental noise is accounted for. They can thus 
be discriminated between on the basis of the genus amplitude. For instance, at the 
0.3° smoothing scale each model is distinguished from the other two at greater than 6-cr 
confidence (Table 1 and we have added all error bars in quadrature). (We again emphasize 
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that the other cosmological parameters, such as h and Qb, are not allowed to vary in these 
models.) These results may be directly compared to the genus measurements from MAP, 
as an independent check of the results inferred from the power spectrum and correlation 
function analyses. 

4.4. Effects of Noise on the Genus 

In order to understand more completely the effects of receiver noise on the genus, we 
can look in detail at what is happening, structure by structure, when noise is added. At a 
high temperature cut, for instance, noise can add a hot spot by adding a small amount of 
positive temperature to a region near the cut; conversely, noise can remove a small hot spot 
by cooling a high-temperature region to below the cut. We have illustrated schematically in 
Figure 7 the possible effects of noise on the genus. In Figure 7b and 7e the noise has neither 
added nor subtracted hot spots, leaving the genus unaltered. In Figure 7a, however, the 
noise has removed a hot spot and reduced the genus by one; in 7c the noise has split a hot 
spot into two and added one to the genus. In Figure 7d the noise has added an artificial 
hot spot and one to the genus while in 7/ the noise has merged two hot spots into one and 
reduced the genus by one. These splittings and mergings can occur in higher multiples. 

We now examine our simulated maps for these effects. In Figure 8a we have plotted the 
u — 1 threshold mask for a subsection of the noiseless OCDM 1.0° FWHM smoothed map. 
In Figure 86 we have plotted the same mask, but for a map with noise added. Several of the 
structures in Figure 8a become subdivided in Figure 86 — noise has split hot spots. Also, 
several new hot spots appear, while at least one cold spot (in the upper right) develops 
within a hot spot. To understand these changes quantitatively, we have counted the number 
of times each case in Figure 7 arises for the OCDM simulations of Figure 8. Results from 
this analysis are given in Table 2. 

Table 2 shows that about 80% (660 out of [78+660+64+15+2+2]) of the structures 
counted in the noisc-frcc map correspond to exactly one structure in the map with noise. 
78 real hot spots occur with no corresponding hot spots in the map with noise, while 
83 (=64+15+2+2) overlap two or more hot spots in the map with noise. Very similar 
correspondence exists in comparing the noise-added map hot spots to the noise-free map 
hot spots, except in the "0" column. About a quarter (249 out of [249+683+13+2+1]) of 
the hot spots in the map with noise had no corresponding hot spots in the noise-free map. 
This is to be expected because the noise adds artificial small-scale bumps on top of the 
signal. These bumps can push a pixel that was just below the threshold (here z/ = 1) to a 
value over the threshold, so that extra hot spots appear (similarly, extra cold spots appear 
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at 1/ = -1). At 1.0° FWHM smoothing about two-thirds (683 out of 1011) of the structures 
in the simulated MAP maps are due to real signal. Therefore most 1.0° scale structures in 
the MAP maps will correspond to real fluctuations in the microwave background. 

5. Conclusions 

In this paper we have generated mock CMB anisotropy maps of representative 
DMR-normahzed CDM cosmogonies appropriate for the MAP satellite experiment. We 
have studied the sensitivity of the simulated MAP data to cosmology, sky coverage, and 
instrumental noise. An accurate knowledge of instrumental noise is essential if the statistical 
tests considered here are to live up to expectations. 

We have focused on the correlation function and the genus statistic as tests of 
cosmogonical models. If the underlying theory has Gaussian random phases, then these two 
statistics are just integrals of the power spectrum. However, since they are computed from 
the data in a very different manner than the power spectrum, they will be useful checks on 
power spectrum computations even in Gaussian random phase models. The genus statistic 
will also be a powerful test of the random phase hypothesis. 

The correlation function statistic provides an accurate measurement of the acoustic 
Hubble radius at decoupling from the simulated MAP data. The Qq = 0.4 OCDM model 
considered here can be distinguished from the spatially-fiat FGDM and AGDM models 
with more than 99% confidence from the location of the acoustic valley in the correlation 
function, even when the expected amount of instrumental noise is present. The FGDM and 
AGDM models can also be distinguished at separations ~ 0.5° by the amplitude of the 
correlation function. Note however that for these models all other cosmological parameters 
(such as Qb and h) are fixed. 

The genus of the GMB fiuctuations can also be accurately measured from the MAP 
data. The genus curves of Gaussian cosmological models allow a 2-a statistical horizontal 
shift of the zero-crossing point by only \Au\ ^ 0.01 (0.04) when the total effective 
smoothing is 0.3° (1°) FWHM. The asymmetry of the genus curve at the positive and 
negative threshold levels is only \Ag/g{u = ±1)| <^ 0.8% (4%) in the Gaussian cosmogonies 
considered. Deviations of the observed MAP data that exceed these values will be evidence 
for non-Gaussian behavior. 

The amplitude of the genus curve, which is a measure of the shape of the power 
spectrum at the smoothing scale, can also be a powerful discriminator between cosmological 
models. Even with the expected amount of instrumental noise and partial sky coverage. 
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the MAP data should discriminate between the cosmogonies considered here at confidence 
levels exceeding 99% just on the basis of the genus amphtude analysis. 

Of possible outcomes from the MAP experiment, a more interesting one would be if 
the results were in close agreement with the Iow-Qq open-bubble inflation model. Not only 
would this tell us the value of but it would affirm infiation and support the open-bubble 
infiation model in particular, thus providing evidence that there are universes other than 
our own (other bubble universes) , all arising out of an initial metastable inflation state. 
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Table 1: Numerical values for genus amplitude per steradian ai u — 1^ 



Smoothing scale 


FCDM 


OCDM 


ACDM 


0.3° 


3650 ± 16 ± 23 


3995 ± 21 ± 40 


3313 ± 14 ±22 


0.5° 


637.6 ± 6.0 ± 13.6 


719.1 ±8.2 ±22.1 


617.9 ± 5.2 ± 14.0 


1.0° 


137.7 ±2.5 ±6.2 


120.9 ±2.6 ±6.6 


149.7 ±2.5 ±6.9 



"^Computed at various smoothing scales for the FCDM, OCDM and ACDM models, from simulated maps 
that account for the expected instrumental noise and the Galactic-plane cut. The first uncertainty listed is 
the statistical error of the fit (1-cr, derived from comparing the genus in eight octants of the sphere), and 
the second uncertainty is that due to cosmic variance. 



Table 2: Effects of noise on the genus^ 





intersecting hot spots 




mapi — > map2 


1 2 3 4 5 


6 


noise-free — > noise 


78 660 64 15 2 2 





noise — > noise-free 


249 683 63 13 2 1 





change of genus 


-1 12 3 4 


5 



'^AIso see Figures 7 and 8. At the threshold v = \, for a given hot spot in either the noise- added map or 
noise-free map, we have counted the number of hot spots in the other map, which that hot spot intersects. 
In the first row, we have the number of times a hot spot in the noise-free map intersects a hot spot in the 
map with noise: noise c;an cause a hot spot to disappear (col. 0; Figure 7a): noise might leave a hot spot 
as is (col. 1; Figure 76); noise might split a hot spot in twain (col. 2; Figure 7c), etc. In the second row, 
we have the number of times a hot spot in the map with noise intersects a hot spot in the noise-free map: 
the noise might have produced a whole new hot spot (col. 0; Figure Id); noise could have left one spot as 
one spot (col. 1; Figure 7e); noise could have merged two hot spots into one (col. 2; Figure 7/), etc. For 
example, 78 hot spots in the noise-free map are eliminated in the map with noise, while 660 of them are left 
as is, and 64 are split into two spots. Similarly, 249 of the hot spots in the map with noise are eliminated in 
the noise-free map, while 683 of them are left as is, and 63 are split into two. The change of genus associated 
with these effects is given in the last row. 
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Fig. 1. — Curves show the CMB anisotropy angular power spectra, as a function of multipole 
of the four models considered: FCDM (sohd), ACDM (dotted), OCDM (dashed), and 
RCDM (long dashed). Tq is the mean CMB temperature. The clouds of dots, tracing the 
curves at small £, are from actual realizations of the a™'s with 0.21° FWHM beam smoothing 
as appropriate for MAP, and show the cosmic variance at each i. The clouds of dots in panel 
(a) do not account for instrumental noise, while those in panels (6) and (c) do. The dots in 
the last two panels correspond to 0.3° and 1.0° total smoothing of the signal, respectively. 



Fig. 2. — Upper panel is a MoUweide projection of the CMB anisotropy map simulated for 
the OCDM model. The map has been convolved by a Gaussian filter with 0.3° FWHM. Lower 
panels are magnifications of the MoUweide projection of the simulated CMB anisotropy maps 
for the OCDM and ACDM models. In each lower panel map the equator runs horizontally 
through the center, where the plots cover 30° horizontally and 24.4° vertically. 



Fig. 3. — Upper panel is the OCDM CMB anisotropy map, in stereographic projection and 
with 1° FWHM Gaussian smoothing, showing the north and south galactic caps {\b\ > 30°). 
The projection is conformal and the shape of structures arc locally undistortcd. Lower 
panel is a similar stereographic projection for the ACDM model. For the purpose of this 
comparison both models have identical phases, in both signal and noise (picked from the 
same random number seeds). 
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Fig. 4. — Correlation functions of the four cosmogonies. Thick curves are the correlation 
functions transformed from the power spectra shown with curves in Figure 1. Thin solid 
curves, tracing the thick ones at large angular separations, account for the convolution with a 
0.21° Gaussian beam in panel (a), with an additional 0.214° Gaussian smoothing in panel (6). 
Dotted curves are the correlation functions from seven independent reahzations of FCDM 
CMB anisotropy skies showing the cosmic variance. Thin solid and dotted curves in panel 
(a) ignore instrumental noise, while those in panel (b) account for it. The Galactic-plane cut 
is not accounted for. Clearly, inclusion of instrumental noise smoothed over a small scale 
(~ beam size) does not change the correlation function much, indicating that the statistical 
uncertainty is dominated by cosmic variance. 
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Fig. 5. — Two-dimensional genus topology from the MAP simulations accounting for 
instrumental noise, with 1.0° FWHM smoothing and for \b\ > 30°. Open symbols are for 
OCDM, closed for FCDM. Error bars are derived from the variance of the genus computed 
independently in each of the eight octants of the sphere. The curves show the functional 
form expected for a random-phase Gaussian field, Nve"^^^'^^ and have been fit to the points 
by adjusting N. The cosmogonies can clearly be distinguished from each other. 
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Fig. 6. — Fit amplitudes (at u — 1) of genus curves for various cosmogonies, as in Figure 
5. Filled squares are results from noise-free maps for FCDM, filled triangles for ACDM, and 
filled circles for OCDM. Open squares are results from maps with noise for FCDM, open 
triangles for ACDM, and open circles for OCDM. The error bars shown in the insets are the 
quadrature sum of the statistical uncertainty and that due to cosmic variance (see Table 1). 
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Fig. 7. — Overlays of simulated maps with and without instrumental noise. This schematic 
illustrates the effect of noise on the genus by comparing how well a given structure in the 
map with noise corresponds to a structure in the noise-free map. If a structure in the map 
with noise, such as a hot spot, intersects exactly one hot spot in the noise-free map, as in 
cases b) and e), then the noise has not altered the genus. However, if the noise is significant, 
artificial hot spots can appear or real hot spots can disappear, both of which would alter 
the genus, a) is a case where the noise has removed a hot spot (there is a hot spot in the 
noise-free map but not in the map with noise), c) where the noise has split a hot spot into 
two. d) is a case where the noise has added a hot spot (there is a hot spot in the map with 
noise but not in the noise-free map), /) where the noise has merged two hot spots into one. 
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Fig. 8 — Threshold masks from a subset (approximately 30° x 30°) of the OCDM CMB 
anisotropy simulations at 1.0° FWHM smoothing: (a) ignoring instrumental noise, {b) 
accounting for instrumental noise (rms noise = 35 /^K per 0.3° x 0.3° area). 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/97 1 1 057v 1 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/97 1 1 057v 1 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/97 1 1 057v 1 



